Quantum Corrections to Dilute Bose Liquids 
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It was recently shown [A. Bulgac. Phys. Rev. Lett. 89, 
050402 (2002)] that an entirely new class of quantum liquids 
with widely tunable properties could be manufactured from 
bosons, fermions and their mixtures by controlling their in- 
teraction properties by means of a Feshbach resonance. Since 
the quantum fluctuations in this regime could in principle 
destabilize these objects, we extend the previous mean-field 
analysis of these quantum liquids by computing the lowest or- 
der quantum corrections to the ground state energy and the 
depletion of the Bose-Einstein condensate and by estimating 
higher order corrections as well. We show that the quantum 
corrections are relatively small, do not lead to the destabi- 
lization of the droplets and are controlled by the diluteness 
parameter n\a\^ <C 1, even though strictly speaking in this 
case there is no low density expansion. 

PACS numbers: 03.75.Hh, 03.75.Nt, 05.30.Jp 



Among the quantum fluids observed so far only '^He 
and *He and atomic nuclei are self-bound/liquid systems. 
R is very likely that this relatively short list could become 
extremely long by the addition of boson droplets (bose- 
lets) , fermion droplets (fermilets) and their mixtures (fer- 
bolets) [1] , which owe their existence to the possibility of 
tuning two-atom interactions and to some unusual prop- 
erties of the three-body system. 

R was predicted about ten years ago [2] and was con- 
firmed experimentally about five years ago [3] that by 
immersing atoms in a magnetic field (and, in principle, 
in electric or laser fields as well) , one can alter essentially 
at will the scattering length a between two atoms. Un- 
der such circumstances, it is possible to create a situation 
where the scattering length is negative and significantly 
larger than the interaction radius. In the regime where 
the magnitude of the scattering length is much larger 
than the typical scale of atomic interactions rg , the three- 
atom system exhibits some very unusual and interesting 
phenomena [4,5]. The three-body observables depend on 
a new scale aa which is not determined by any two-body 
observable. In the language of effective field theory, this 
corresponds to the fact that a three-body force appears 
at leading order of the low energy expansion [6]. Also, 
three-body observables are periodic functions of ln(a/a3) 
and the three-particle scattering amplitude can be made 
arbitrarily large by fine-tuning a. By choosing a neg- 
ative, large (|a| 3> tq) and with a three-body bound 
state close to the threshold, a rather unique situation is 
created: two low energy atoms experience an effective 



attraction, but at the same time three atoms will expe- 
rience an effective repulsion. This situation is very well 
understood theoretically and it is known as the Efimov 
effect [4,5]. A large ensemble of such atoms can condense 
in the same way as water droplets condense from vapor. 
Since such condensation can occur for any atom species 
for which a Feshbach tunable resonance exist, an entirely 
new class of (bose, fermi, or mixtures) quantum liquids 
[7] can be created. Their equilibrium density is deter- 
mined by the interplay between the two-body attraction 
and the three-body repulsion between particles in the 
Efimov regime described above. Since the magnitudes of 
the two-body attraction and of the three-body repulsion 
are to a significant extent under experimental control, the 
basic properties of these new quantum liquids are widely 
tunable [8]. 

In Ref. [1], the basic properties of this new class of 
quantum fluids were established in the mean field ap- 
proximation. It is not at all obvious that quantum fluc- 
tuations could not destabilize a boselet. It is generally 
expected that the mean field approximation is rather 
accurate for dilute systems [9-11]. Under normal cir- 
cumstances (when a > and a ~ tq) it is sufficient 
to take into account only the binary atom-atom colli- 
sions and the effects of quantum fluctuations, whereas 
the triple atom and higher collisions are rather small and 
controlled by the diluteness parameter \/no? <^ 1 [9-11]. 
The regime a < was considered until recently as in- 
trinsically unstable towards collapse, unless the number 
of bosons is smaller than about 1,500 or so in a trap 
[9,12]. In the specific Efimov regime, we are interested in 
the fact that the two-body scattering amplitude is large 
(a < 0, \a\ ^ ro), but the three-body amplitude is even 
larger. Even though the system is dilute with respect 
to the two-body collisions (n|a|'^ ^ 1), it is not obvious 
that diluteness with respect to the three-body collisions 
is achieved and one might suspect that the mean-field 
approximation could be violated. It is imperative to un- 
derstand the character and the magnitude of these correc- 
tions to the ground state properties in the case of these 
new quantum liquids and determine whether quantum 
fluctuations can lead to instabilities. One can naturally 
expect that such effects would perhaps manifest them- 
selves particularly strongly in the case of bosons and less 
so in the case of fermions. For this reason, we shall focus 
our attention here on Bose systems only. 

We start with the Lagrangian density 
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where tp^ and are the creation and annihilation oper- 
ators for a boson, jj is the chemical potential, and the 
coupling constants are 
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Here a is the two-body scattering length and sq ~ 
1.00624. dx and ^2 < are universal constants, whose 
numerical values have been determined recently [13]. The 
values of the coupling constant are determined by con- 
sidering the scattering of two and three particles at zero 
momentum. 03 is the value of the two-body scatter- 
ing length for which a three-body bound state has ex- 
actly zero energy. Unlike di^2 and so, the parameter 03 
is system dependent and is also a genuine three-body 
characteristic. The rest of the symbols have their usual 
meaning. We have introduced a new dimensionless quan- 
tity T » 1 in Eq. (3), which will prove very conve- 
nient in the power counting. We will only use the La- 
grangian in Eq. (1) for particle momenta much smaller 
than %/\a\ so that it is legitimate to subsume the com- 
plicated dynamics on the scale fi/a, leading to the Efi- 
mov effect into a contact three-body interaction. We 
show below that the typical loop momenta are of order 
Q ~ 7i/(|a|T^/^) <C ?i/|a|. More precision can be sys- 
tematically attained by including terms in Eq. (1) with 
more derivatives or fields, but their effect is strongly sup- 
pressed. As wc mentioned above, the interesting regime 
is when 52 < and 53 > [1]. 

Since we are interested in the condensed state of a 
system of bosons, it is natural to split the field ip into a 
classical (c-number) condensate and the fluctuations ip: 
tp ^ ip + (f). {(f) can be chosen real in this case.) In terms 
of the new variables, the Lagrangian density becomes 
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where x = Qi'i?' +530^- The interaction terms include 
vertices with three ^''s (~ g2<P and ~ 530"^), four i/i's (~ 
and ^ 530^), five ■0's (~ 93<^-, and six ^''s (~ 33) and nor- 
mal ordering is implied as well. At tree level/mean field, 



the particle number density is 



and the energy 



density and chemical potential (determined by imposing 
the condition = 0) are given by 
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In writing down the contribution quadratic in i/)'s in Eq. 
(4) wc have implicitly taken into account the condition 
(6). At equilibrium/zero pressure {P = nd£/dn — £ = 0) 
the values for the particle number density, energy density, 
and chemical potential are easily determined to be 
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Notice that only at equilibrium the chemical potential 
becomes equal to the energy per particle, see Eq. (9). 
The tree level/mean field approximation is valid when 
the system is dilute, which holds as long as T S> 1. At 
and near equilibrium, the contribution of the two-body 
and three-body collisions are comparable in magnitude. 
The three body collision term, which in dilute systems 
is typically negligible, is important in this case, because 
the three particle bound state is close to the threshold 
and the particles arc in the Efimov regime. 

From the Lagrangian (4) we can determine the normal 
and anomalous propagators for the Bogoliubov quasipar- 
ticle, which can be conveniently represented as a matrix 
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where cuk is the Bogoliubov dispersion relation for quasi- 
particle excitations. 
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Since 32 < 0, the sound velocity s = \fxirn for long 
wavelengths is imaginary for small particle number den- 
sities, namely when 
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and the homogeneous matter at these densities (n < n^) 
is unstable towards collapse and condensation into denser 
droplets. For densities Ug < n < no, the pressure is neg- 
ative and such a system tends to increase its density and 
binding energy. If the density is n > no, the internal 



2 



pressure is positive and the system tends to expand to- 
wards equlibrium and lowers its energy, unless external 
walls exert inward pressure. Unlike the phonons in liquid 
^Hc, the curvature of the quasiparticlc excitation LOk as 
a function of the wave vector k is positive. As a result 
these quasiparticles have a finite lifetime due to decay 
into two quasiparticles, corresponding to a width T ^ 
[11,16], since the velocity of a quasiparticlc with finite 
wave vector k is larger than the sound velocity. 

The quantum fluctuations can now be evaluated at the 
one-loop level with usual techniques. We have chosen to 
use dimensional regularization for the evaluation of var- 
ious diverging integrals [17]. Old fashioned techniques 
[9-11], such as diagonalization of the quadratic part of 
the Lagrangian (4) with the a Bogoliubov transforma- 
tion and subsequent evaluation of the corrections to the 
ground state density and energy would lead to identical 
results. The particle density and the energy density at 
one-loop level are given by the following expressions 
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The first correction term to the energy density (and sim- 
ilarly to the number density as well) is formally identical 
to the Lee and Yang term [9-11], which was computed, 
however, for a Bose gas with repulsive two-body inter- 
action. The equilibrium values for the particle number 
density (both condensed and noncondcnsed), the energy 
density, and energy per particle, after including first- 
order quantum corrections become 

1 1 ^ 72 \ ^ 23/2 



(16) 



Surprisingly, Eqs. (13) and (14) determine correctly both 
the linear and the quadratic corrections in a for energy 
per particle £/n. Therefore, the first order quantum cor- 
rections arc controlled by the same diluteness parameter 
^/n\a^ = T~^/2 ^ 1 as in the case of a Bose gas with 
repulsive interaction. In spite of this formal apparent 
similarity, the meaning of the present result for the con- 
densed Bose liquid when a < is qualitatively different. 
There is strictly no low-density expansion in the present 
case, since for densities n < Us = 2no/3, the system 
is unstable towards long-wavelength density fluctuations 
as the sound velocity is imaginary, see discussion of Eqs. 
(11) and (12). 

A general argument can be given for the suppression of 
higher loop corrections. Consider a generic diagram con- 
tributing to the energy density (that is, without external 



legs) containing I propagators, L loops and rii vertices 
with i legs. It can be estimated by 
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where Q ~ y^mx ~ l/dc^lT^/^) is the typical loop mo- 
mentum [18]. When using dimensional regularization no 
powers of the cutoff are present and the only momentum 
scale in these diagrams is set by ^Jvvx- Combining the 
above estimate with the relations 
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we find 



(Energy Diagram) ~ m^x*"^^^^"'^ ~ E^T'^I'^. (20) 

Consequently, up to logarithmic corrections, all higher- 
order loop corrections are suppressed in the dilute limit 
by powers of \/n|aP = T"^/^ ^ 2. As in the case of 
dilute Bose gases with repulsive interactions [9-11,14,15] 
this parameter is proportional to the ground state frac- 
tion of non-EEC particles. In this respect, this new class 
of quantum liquids, in particular the boselets, is quali- 
tatively different from liquid ^He, where the fraction of 
non-BEC particles is close to 90%. 

For large T 1, the self-bound system is dilute. The 
states with positive pressure, however, can be dense. 
These states could still be dilute with respect to two- 
body collisions, however of such a density that |g2| <C 
(7371. In this regime, the loop expansion breaks down 
and a further rcsummation is required. The energy per 
particle of the system is large and the system can ex- 
ist only if external pressure is applic^d. It is very likely 
that the system undergoes a transition to a coherent 
mixture of atoms and trimers, similar to the situation 
put in evidence in gaseous BEC, where the coexistence 
of an atomic BEC with a molecular BEC was predicted 
[21] and experimentally observed [22]. The formation of 
trimers will have two effects, lowering the particle den- 
sity and lowering the energy. The spatial size of a trimer 
when a < is smaller than |a|, even if the trimc;r bind- 
ing energy is vanishing. The wave function of an Efimov 
state of zero energy is normalizable and the wave function 
is concentrated predominantly at interparticle distances 
\vi — Vj\ < |a|, where = 1,2,3 are particle labels in a 
trimer [19]. Since the average separation between atoms 
is presumed to be larger than the scattering length, when 
all atoms convert into trimers, the trimer number den- 
sity decreases roughly by a factor of 3 when compared 
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to the atomic number density [1]. The trimer phase is 
most hkely stable, provided the trimer-trimer scattering 
length is positive, which is most probable [9,20]. When 
a trimer is formed, a significant amount of energy is re- 
leased as well. These two factors, the significant drop 
in density and energy, should be explicitly taken into ac- 
count when evaluating the energy density at higher parti- 
cle densities, where triple collisions might dominate. The 
self-bound state we have considered is only meta-stable. 
It can decay either by a four particle process, leading 
to the formation of a trimer and one energetic atom, or 
through a three-body process with the formation of a 
deeper dimer and one energetic atom (the formation of a 
shallow dimer is excluded since a < 0). The rate for the 
first process is proportional to %n^aJ /m and likely very 
small. The formation of a deep dimer was considered in 
Ref. [23] and it was found to be of the order of Tin^a^ /m. 
This decay rate is, however, enhanced when a three-body 
state is close to the threshold, where the correct rate can 
be determined only through a non-perturbative calcula- 
tion, but this is not yet attempted. This recombination 
rate is related to the widths of the Efimov states due to 
the existence of deep two-body bound states, which fol- 
low the same exponential behavior as the energies of the 
Efimov states [24]. 

We thank G.F. Bertsch and D.B. Kaplan for discus- 
sions. 

Note added. Boselets made of spin polarized tritium 
atoms are the likely candidates, since the regime we study 
here, when 52 < and 33 > 0, could apparently be easily 
reached according to Ref. [25]. In this particular case, 
three-body recombination processes are absent. 

Note added in proof. Recently Braaten and Ham- 
mer [26] have performed a nonperturbative calculation 
of three-body recommbination rates in dilute gases. 
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